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Presenting a general phase approach to stochastic processes we analyze in particular the Fokker- 
Planck equation for the noisy Burgers equation and discuss the time dependent and stationary 
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The strong coupling aspects of systems driven stochas- 
tically far from equilibrium present a formidable chal- 
lenge in modern statistical physics and soft condensed 
matter. The phenomena in question are ubiquitous and 
include turbulence in fluids, interface and growth prob- 
lems, chemical reactions, self organized critical systems, 
and even economical and sociological models. 

In recent years much of the focus of modern statis- 
tical physics and soft condensed matter has shifted to- 
wards such systems. Drawing on the case of static and 
dynamic critical phenomena in and close to equilibrium, 
where scaling, critical exponents, and universality have 
served to organize our understanding and to provide cal- 
culational tools, a similar approach has been advanced 
towards noncquilibrium phenomena with the purpose of 
elucidating scaling properties and more generally the 
morphology or pattern formation in a driven state. 

Whereas perturbative field theory together with the 
dynamic renormalization group have proven successful in 
the context of dynamic critical phenomena, the extension 
to nonequilibrium phenomena is often plagued with both 
technical and conceptual problems. This is related to the 
occurrence of strong coupling features encountered for ex- 
ample in the notable case of hydrodynamical turbulence, 
and there is a need for the development of appropriate 
nonperturbative theoretical tools in order to access the 
strong coupling regime. 

In this context the noisy Burgers equation for the slope 
u n — W n h (n = 1, ..d) of a growing interface jjj, 



du n 
~dt 



XUpVpUn 



(1) 



or, equivalcntly, the Kardar-Parisi-Zhang (KPZ) equa- 
tion M for the height h, 



dh „ 9 , X„ ,„ , 
— = vV 2 h + -V n hW n h + rj 
at 2 



(2) 



provide maybe the simplest continuum description of an 
open driven nonlinear system exhibiting strong coupling 



features such as dynamical scaling and pattern formation. 
In ([!]) and (|^) v is a damping constant or viscosity char- 
acterizing the linear diffusive term, A a coupling strength 
for the nonlinear mode coupling or growth term, and r\ a 
Gaussian white noise driving the system into a stationary 
state and correlated according to 



(r)(x n ,t)r](x' n ,t')) = A JJ(5(x„ - x' n )S(t - t') 



(3) 



characterized by the noise strength A. 

Notwithstanding the simple form of ([!]) and (^), the 
driven Burgers equation introduced originally in order to 
model aspects of turbulence and the KPZ equation pro- 
viding the simplest description of a growing interface, the 
morphology and scaling properties embodied in (|l]) and 
(Q) have been difficult to extract and a full understand- 
ing of (0) and (|) remains one of the important issues in 
nonequilibrium statistical physics || . 

Besides perturbation theory in A ^ which as regards 
the scaling properties provides the roughness and dy- 
namic exponents (Q,z) = (1/2,3/2) in d = 1, but, other- 
wise, is limited to an e expansion about the (lower) crit- 
ical dimension d = 2, yielding a kinetic phase transition 
above d = 2 separating a weak coupling phase (the A = 
universality class) with exponents (£, z) = ((2 — d)/2,2) 
from a strong coupling phase and to all orders in e the 
exponents (£, z) = (0, 2) on the phase line ||, nonpertur- 
bative methods include 1) in the d = 1 case mapping to 
spin models Q and information gained from lattice mod- 
els [Q, 2) mapping to directed polymers in combination 
with replica methods ||, 3) mode coupling expansions 
and, most recently, 4) operator expansions yielding 
the strong coupling exponents ((, z) = (2/5, 8/5) in d = 2 
and (C, z) — (2/7, 12/5) in d = 3, corresponding to skew- 
ness in the height distributio n |H | . 

In a recent series of papers [ [LOf we advanced a nonper- 
turbative approach to the noisy Burgers equation which 
purports to elucidate both the morphology and scaling 
properties of a growing interface in d — 1. Arguing that 
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the noise strength A is the relevant nonperturbative pa- 
rameter, driving the system into a stationary state, and 
thus circumventing 1) the limitations of perturbation the- 
ory which is based on an effective expansion in A 2 A/i/ 3 
assuming regularity in A M and 2) the self consistency 
assumptions underlying the mode coupling approach || , 
the method was based on a weak noise saddle point ap- 
proximation to the Martin-Siggia-Rose functional for- 
mulation |ll]] of the noisy Burgers equation ([[]). Im- 
portantly, the method yields coupled deterministic field 
equations for the slope u and a noise field (p (character- 
izing 77), replacing the stochastic Burgers equation, ad- 
mitting soliton solutions and as a result a many body 
formulation of the pattern formation of a growing inter- 
face in terms of a dilute gas of propagating solitons with 
superposed linear diffusive modes. The canonical form of 
the approach also yields 1) the gapless soliton dispersion, 
E oc Xp z , z = 3/2, and gapful diffusive mode dispersion 
u> = —iv{k 2 + (Xu/2iy) 2 ), u is the soliton amplitude, 2) 
recovers the scaling exponents (£, z) = (1/2,3/2) and an 
expression for the scaling function, and 3) associates the 
Burgers universality class with the leading gapless soliton 
excitation. 

In the present letter we develop a general canonical 
phase space approach to a stochastic Langevin equa- 
tion of the Burgers type with additive white noise. This 
method which emerged from our previous studies alluded 
to above allows us to discuss and in some cases derive the 
stationary and time-dependent weak noise solutions of 
the associated Fokker-Planck equation for the probabil- 
ity distributions. In particular for the Burgers equation 
(|l|) the time-dependent and stationary distributions are 
given by 



P{u n , T) oc exp 



A 



S(u n ,T) 



P st (u n ) oc lim exp 

T — >oo 



A 



S{u n ,T) 



where the action has the canonical (symplectic) form 



S = £ d d xdt(p n ^-H 



(4) 
(5) 

1 

(6) 



with Hamiltonian density 

U = p n {vV 2 u n + \u m V m u n - (l/2)V„V m p m ) , (7) 
yielding the coupled Hamiltonian equations of motion 
d 



dt 
d_ 
dt 



- Au m V m u„= vV u n - V„V m p r , 



(8) 



- Aw m V m ) Pn= Pn 



u m ) . (9) 

The above equations allow an analysis of the time- 
dependent distribution P(u n ,T). In principle we have 



to solve the canonical field equations (||) and (Q), deter- 
mining the orbits in p n u n phase space, and compute the 
action (||) and thus P{u n ,T) according to (|j). 

The general framework is developed along the follow- 
ing lines: The Fokker-Planck equation pertaining to a 
general Langevin equation with additive noise for the 
stochastic variable qi (i is a discrete and/or continuous 
index) , 



-dT = -2 Fn{qi) + 7]n ' 

(Vn(t)Vm(t')) = AK nm S(t-t') , 

has the form, denoting V n = d/dq n , 
OP 1 

— = - V n [AK nm V m P + F n P] 
Searching for a solution of the form 



P oc exp 



1 



(10) 
(11) 

(12) 
(13) 



it is an easy task to show that to leading order in the noise 
strength A, the action S(q n ,t) satisfies the Hamilton- 
Jacobi equation p2[ 



r99 

— + H(q n ,V n S) = 



(14) 



where the energy E — H and the canonically conjugate 
momentum p n — V n S. The Hamiltonian H has the gen- 
eral form 



H = (l/2)(K nm p n p m - F n p n ) , 
implying the Hamiltonian equations of motion 

^-K v --F 



dp n 

dt 



1 



(15) 

(16) 
(17) 



Assuming for simplicity that F n — » for q n — ► the en- 
ergy surfaces have the characteristic submanifold struc- 
ture depicted in Fig. 1. 

The origo in phase space constitutes a hyperbolic sta- 
tionary point defined by the unstable zero-energy sub- 
manifold p n = 0, the transient manifold, and, assuming 
the existence of a stationary state, a stable submanifold 
defined by K nm p m — F n orthogonal to p n , the stationary 
manifold. The stationary state is determined by orbits on 
the zero-energy manifolds whose structure thus charac- 
terizes the nature of the stochastic problem. The action 
S and hence the distribution P are given by 



S(q n ,T,q' n 



dt 



dqn 
'■ dt 



H 



P(qn,T,q' n )(x exp [S(q n , T, q' n )/A} 



(18) 
(19) 



where the orbit from q' n to q n is traversed in time T. 
The stationary distribution P s t{q n ) is thus obtained in 
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the limit T — > oo and E = 0, assuming that E(T) oc 
exp [—const.T] for T — > oo, 



ist(ffn) = exp 



1 /°° ,. dq n 



(20) 



i.e., P s t is determined by the infinite time orbits on the 
zero- energy manifold. 




FIG. 1. Canonical phase space in the general case. The 
solid curves indicate the zero-energy transient submanifold (I) 
and stationary submanifold (II). The stationary saddle point 
is at the origin. The finite time (T) orbit from q' n to q n mi- 
grates to the zero-energy submanifold for T — > oo. 

This structure of phase space permits a simple non- 
stochastic, deterministic discussion of the approach to 
the stationary state of a damped noise-driven system 
in terms of dynamical system theory. Referring to Fig. 
1, consider an orbit from q' n to q n on the energy sur- 
face E(T) traversed in time T. In order to attain the 
stationary state E(T) — > in the limit T — > oo. For 
E ~ the initial part of the orbit moves close to the 
p n = submanifold and from (|l(]) is determined by 
dq n /dt = —(l/2)F n i.e., the deterministic noiseless ver- 
sion of the Langevin equation (fL0|). In the absence of 
noise the motion is transient and damped. The orbit 
slows down near the stationary point (the origo in phase 
space) before it picks up again and moves close to the 
other submanifold K nm p m — F n _L p n . This final part of 
the orbit terminating in q n at time T thus corresponds 
to the establishment of the stationary state. Markovian 
behavior, i.e., the independence of the initial configura- 
tion, is associated with the long (infinite) waiting time 
near (at) the stationary point. 

In the special case K nm — S nm and F n = V„$, corre- 
sponding to an effective fluctuation-dissipation theorem 
and an underlying free energy <I>, the energy ([l5|) and 
equations of motion ( |l6| ) and ( |l7| ) are consistent with the 
zero-energy submanifolds p n — and F n — p n and yield 
the equilibrium distribution P s t(ln) exp [— $/A]. In 
the general case a determination of the time-dependent 
and stationary distributions require a knowledge of the 
energy submanifolds in combination with a solution of 
the canonical equations in order to determine the orbits. 

The noisy Burgers and KPZ equations (Q) and ^ fall 
within the scope of the general framework summarized 



above. With the identification qi (t) — > u n (x m ,t), K nm — > 

^ 2 Tl m H x m and F n -> -2(^V 2 U„ + Au m V m M„), 

we obtain (^)-(g). Note that the canonical momentum, 
the 'noise field' p n , is essentially a 'slaved variable'. 

In the case d = 1 which is the basis for our discussion 
here, the equations (H) and (||) reduce to the form 



- AuV u= v\I 2 u - V 2 p , 



0_ 

dt 



- XuV ) p= -vV 2 p 



(21) 
(22) 



The distribution (jij) is determined by the action 
which has the general structure 

S = S st {u) + S si , cw (u,T) + S so i{u,T) . 



(23) 



The stationary distribution P s t (u) given by S s t{u) is 
easily found by noting that ( ^l|) and (|22j ) coincide on the 
zero-energy submanifold p = 2vu (the energy density Ti. 
becomes a total derivative yielding E = 0). Inserting 
p = 2vu and E = in @ and (g) in the limit T -> oo 
and integrating over time we obtain the the symmetric 
Gaussian stationary distribution |13|| 



P si (u) oc exp 



(v/A) / dx u(xf 



(24) 



Also, setting p = 2v(u+5u) we find to leading order in 8u 
(d/dt- XuV)Su = uV 2 8u. Noting that d/dt-\uV is in- 
variant under the Galilean transformation: x — * x — Xuot, 
u — > u + uq, and choosing an instantaneous frame with 
vanishing u, Su oc exp [— vk 2 t], implying that the or- 
bits approach the zero-energy stationary submanifold 
p = 2vu. 

The long time skew distribution Ps^^^u^T) is de- 
termined by S' s kew(M,T). For A = we obtain in 
wave number space the symmetric contribution, Uk = 
J dx exp(—ikx)u(x) , 



aO 

skew 



,{u k ,T) =-v I —\u k \ 2 exp[-2^ 2 r] , (25) 



defining a crossover time T co oc 1/vk 2 . For a finite sys- 
tem k oc 1/L (L is the system size), i.e., T co oc L 2 /v, 
yielding the dynamic exponent z — 2 in accordance with 
the diffusive mode contribution. 

For A 7^ and for large T approximating the orbit 
close to the manifold by an orbit on the manifold, in- 
serting p = 2vu in (^l|), we obtain the deterministic 
Burgers equation with rj = and viscosity —v which 
can be solved by means of the Cole-Hopf transformation 
Q. Thus setting u = V/i and h = — (2i// X)hiw yields 
the diffusion equation dw/dt — —v~S/ 2 w for w which 
is solved by means of the Green's function G X (T) = 
[47r^T]- 1 / 2 exp[-a; 2 /4i/r]. We obtain S skew (u,T) = 
v J dx u'(x) 2 , where u — Vh and u' = V/i' are related 
according to the nonlinear expression 



3 



exp[-(A/2z/)/i'] = / dx'G x - m ,{T)exv[-(\/2v)h] , (26) 



P sol (/i,T) oc exp[-(^/A)(l/AT) 1 / 2 / l 3 / 2 



(33) 



giving rise to the distribution 



P(u,T) oc P st (w)exp 



dx u'(xY 



(27) 



which by inspection is skew. To order A in wave number 
space we also have, setting G k ,r — exp[— z/fc 2 T], 



P{u k ,T) cx P st (u k )P^ cw (u k ,T)P s i cw (u k ,T) 



(28) 



where P s ° kcw oc exp[-S s ° kew /A], P s A kow « exph^skew/AL 
and 



2tt 2tt 

Fk,k',T — [G k ,2T — G k ,TG k + k ' ,rG k < ,T]u k U~ k - k >h k i , 

(29) 

exhibiting skewness. 

The short time distribution P so i(u,T) given by 
S so \(u,T), corresponding to an orbit off the zero-energy 
manifold, is determined by the soliton-diffusive mode 
contribution discussed in |10|. For a single soliton with 
boundary values u + and u_ the propagation velocity v 
is given by the soliton condition 



-2v/X 



(30) 



However, only the left hand soliton (u + < u_) carries 
nonvanishing energy, momentum and action according to 
the assignment E = (2/3)^A(w 3 - w 3 ), n = v(u\ - m^), 
and S = (l/6)z^AT|u + — u_| 3 (note the Galilean invari- 
ance of S). The action of a multi-soliton configuration 
constituting a growing interface is then given by 



S so i{u,T) = ^Ai/T^ 

' lhs 



(31) 



where summation is over left hand solitons (lhs) only. 

Owing to the constraint imposed by the soliton con- 
dition (pG) and the non-integrability of the equation of 
motion we can only give a qualitative discussion of P so i- 
Since the saturation width of an interface is a finite size 
effect time scale separation only occurs for a finite sys- 
tem. Noting that the propagation of solitons and the 
imposition of periodic (or bouncing) boundary condi- 
tions, in order to ensure growth in h, endows the ve- 
locity with a scale, i. e., v ~ L/T, we obtain, insert- 
ing it + + U- = —2v/\ ~ u + — w_ in ( |3l"| ) and from 
Pol ~ exp[— Ssoi/A] the soliton crossover time 



T c s o1 oc (1/A)(^/A) 1 / 2 L 3 / 2 



(32) 



which is also is consistent with the dimensionless ar- 
gument A(A/i/) 1 / 2 i/a; 3 / 2 in the scaling function for the 
slope correlations discussed in [p|,^0|. 

In the short time regime T <C T^° l the soliton con- 
figurations contribute to P. Noting that \u + — u_| 3 ~ 
(mL) 3 / 2 (TA)~ 3 / 2 - h 3 ^(TX)-^ 2 we obtain, inserting in 
(|3l"|) (h is measured relative to the mean height) 



in accordance with the directed polymer-replica based 
result H and the exact results for the asymmetric exclu- 
sion model [[uj. The skewness of the distribution then 
arises from the bias in the statistical weight exp[— S/A] 
assigned to the left and right hand solitons giving rise to a 
predominance of right hand solitons (S — 0), correspond- 
ing to relative forward growth. In the long time regime, 
T ^> P c s °', the soliton contribution vanishes and only the 
diffusive modes and their interactions contribute to P. 

In this letter we have outlined a novel canonical phase 
space approach to Langevin equations with additive 
noise; details will appear elsewhere. In addition to pro- 
viding insight from dynamical system theory the method 
also yields a calculational tool for the determination of 
the weak noise probability distributions. In particular 
we have applied the method to the noisy Burgers equa- 
tion in d = 1 and derived expressions for the skew finite 
time probability distribution. In the short time regime 
our heuristic result is in agreement with the directed 
polymer-replica method. 

Discussions with M. Kosterlitz, T. Hwa, J. Hertz, P. 
Cvitanovic, K.B. Lauritsen, and A. Svane are gratefully 
acknowledged. 
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